Quantum diffraction and interference of spatially correlated photon pairs 

and its Fourier-optical analysis 



O 
O 

<N 



<N 

> 

O 

CO 

o 

\o 
o 



Ryosuke Shimizu, 1 Keiichi Edamatsu, 2, 1 and Tadashi Itoh 3 ' 1 

1 CREST, Japan Science and Technology Agency, J^-l-8 Honcho, Kawaguchi 332-0012, Japan 
2 Research Institute of Electrical Communication, Tohoku University, Sendai 980-8577, Japan 
3 Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan 

(Dated: ) 

We present one- and two-photon diffraction and interference experiments involving parametric 
down- converted photon pairs. By controlling the divergence of the pump beam in parametric down- 
conversion, the diffraction-interference pattern produced by an object changes from a quantum 
(perfectly correlated) case to a classical (uncorrelated) one. The observed diffraction and interference 
patterns are accurately reproduced by Fourier-optical analysis taking into account the quantum 
spatial correlation. We show that the relation between the spatial correlation and the object size 
plays a crucial role in the formation of both one- and two-photon diffraction-interference patterns. 
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INTRODUCTION 



Optical engineering is one of the most precise, sophis- 
ticated, and valuable technologies. However, its guid- 
ing principle is still based on "classical wave optics" 
Recently, the non-classical nature of light has attracted 
much attention because of its possible application to 
novel optical technologies, including quantum informa- 
tion processing Q and quantum imaging [3|. Actually, 
over the past few decades, a considerable number of stud- 
ies have been conducted on non-classical light, which can- 
not be explained by classical wave optics. In such studies, 
photon pairs produced by spontaneous parametric down- 
conversion have frequently been used because they can be 
easily generated and can provide a superior environment 
to realize the concept of the non-classical nature of light. 
With the development of studies on entangled photons, it 
has become clear that entanglement, or quantum correla- 
tion, between the constituent photons plays a crucial role 
in the formation of diffraction patterns. Klysyko and his 
co-workers have carried out pioneering theoretical work 
in this field Q. Based on this proposal, Strekalov and his 
co-workers have reported non-classical effects in diffrac- 
tion phenomena, known as "ghost-interference" |5| . This 
work has been developed into one of the fields of quan- 
tum imaging. Meanwhile, the concept of the "photonic 
de Broglie wave" has been proposed as an intuitive and 
essential way to understand the interferometric proper- 
ties of a multiphoton state Q. Within this concept, the 
photonic de Broglie wavelength of an ensemble photons 
with wavelength A and number of photons n can be mea- 
sured to be X/n. Recently, a number of experimental 
demonstrations of this phenomenon have been realized 
EI • This reduced interferometric property 
of a multiphoton state has attracted great interest as a 
light source for novel imaging technology, achieving high- 
resolution imaging at sub-wavelength resolution. This 
proposal is known as "quantum lithography" 0, H^ . 
D'Angelo and her co-workers have shown, as a proof-of- 



principle experiment of quantum lithography, that the 
width of the two-photon diffraction patterns is narrower 
than that of classical light by a factor of 2 For 
its simplicity, Young's double slits have been frequently 
used in these experiments, and the expected counting 
rate was obtained by calculating a quantum-mechanical 
fourth-order correlation function. The use of Young's 
double slit constitutes a simple and thoughtful approach, 
but the method of directly calculating the fourth-order 
correlation function is relatively complicated for analyz- 
ing arbitrary patterns of two-photon diffraction. In this 
paper, we present a simple method using Fourier-optical 
analysis for a two-photon state, considering the quantum 
correlation between the constituent photons. This ap- 
proach corresponds to the Fraunhofer diffraction of the 
classical optics, and is applicable to any arbitrary ob- 
ject. To demonstrate the essence of this concept, we 
measure the diffraction-interference patterns of paramet- 
ric down- converted photons through a transmission grat- 
ing, as an example of arbitrary objects, by means of both 
one- and two-photon detection schemes. We show that 
the observed patterns can be quantitatively explained by 
Fourier-optical analysis. 



II. FOURIER-OPTICAL ANALYSIS OF A 
TWO-PHOTON STATE 

In the process of parametric down-conversion, the 
properties (energy, momentum, and their uncertainty) 
of signal and idler photons are determined by the phase 
matching condition in a nonlinear crystal, concerning the 
length and dispersion of the crystal [l5| and the diver- 
gence of pump beam In the following, we assume 
monochromatic, paraxial, and thin crystal approxima- 
tion for the sake of simplicity. Under these conditions, we 
need to consider only the transverse components of the 
wave vector 0. Thus, hereafter, we discuss diffraction- 
interference in a one-dimensional system. 

Although the actual form of a two-photon state can 
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be derived from phase matching condition of parametric 
down-conversion 0, ^3 ? providing an explanation for 
two-photon diffraction phenomena by an arbitrary two- 
photon source, we begin to consider the following two- 
photon state in relation to an object: 



I dx I 



dx i dx'F(x,x')ai(x)ai(x')\0) 



(1) 



where F(x,x') represents the two-photon probability am- 
plitude at the object, and a\x) is the creation operator 
of the photon at a position x. Since the two photons are 
in the same mode, the two-photon amplitude F at the 
object should have a symmetrical form. Here, we assume 
a spatially correlated form of the two-photon amplitude: 

F(x, x) = i {A(x)B(x') + A(x')B(x)} G(x - a?'), (2) 

where A(x) and B(x') are the probability amplitude 
of photons at the position x and respectively, and 
G(x — x') represents the symmetrical spatial correlation 
between them. For instance, photons from the paramet- 
ric down- conversion process, which is used in our exper- 
iment, naturally form such a state, because a signal and 
its conjugate idler photons generate in almost the same 
position with some uncertainty in a nonlinear crystal. 
The two-photon counting rate at the object is given 
by the second-order correlation function: 



R {2 \x,x') = {^\a ] (x')a ] (x)a(x)a(x') \ 
= \(0\d(x)d(x')\^)\ 2 
= \F(x : x')\\ 



(3) 



where a(x) is an annihilation operator of the photon at 
the position x. In the same manner, the two-photon 
counting rate at the far field can be written as a function 
of the transverse wave number q: 



R (2 Hq,q') = (4>\aHq')aHq)a( q )a(q')\ 
= \(0\a(q)a(q') 



(4) 



where a){q) and a(q) are the creation and annihilation 
operator of the photon with the transverse wave number 
q, respectively. The operator a(q) at the far field is given 
by the Fourier-transformation of the operator a{x) at the 
object: 



1 f ^ 

l(q) = J a(x) exp(iqx)dz 



(5) 



Thus the two-photon amplitude at the far field is given 

by 

F(q,q') = (0\a(q)a(q')\^) 



— J dx J dx (0| d(x)d(x') 
x exp [i(qx + q'x')\ 

— J dx J dx' F(x, x) exp [i(qx + q'x')] . 



In our case, A(x) = B(x) because the two photons pass 
through the same object. Therefore Eq. © is reduced to 



F(q, q') = ±-jdxj dx'A(x)A(x')G(x - x') 

x exp [i(qx + q'x')] . 

This expression is exactly the same as Fraunhofer diffrac- 
tion for a two-photon wave packet except for the intro- 
duction of the function of G(x — x'). Unlike other Fourier- 
optical approaches to two-photon interference 0, H^ . 
we treat the spatial correlation of the photon pairs as a 
phenomenological correlation function. As will be shown 
later, the transverse correlation function G(x — x') plays 
a crucial role in the formation of diffraction-interference 
patterns. The counting rates for two-photon (R^) and 
one-photon (i^^ 1 ^) detection at the object are given by 







R {2 \q,q')= F(q,q>) 



R {1) (q) = J dq'R^(q,q'). 



(8) 
(9) 



First, we consider two extreme cases in which we can 
analytically calculate Eqs. (|HJ) and ©. One is the case in 
which the two photons are emitted without any spatial 
correlation, corresponding to G(x-x') = const., i.e., the 
classical case. The other is the case in which the two 
photons have a perfect spatial correlation, corresponding 
to G(x — x') = S(x — x'). In the former case, there is no 
transverse correlation between the two photons, and the 
Eq. can be rewritten as follows: 



F(q,q') 



h J dxA{x) 



ex.-p(iqx) 



dx' 'A(x') exp(iq'x') 



T[A]{q)F[A]{q'), 



(10) 



where J~[A] represents the Fourier transform of A. In 
this case, the two-photon amplitude is separable into two 
independent photon amplitudes. Therefore, the two- and 
one-photon counting rates are given by 



R^ ( q , q >) oc \F[A]( q )\ 2 \F[A}(q') 
cx R^(q)R^(q'), 



R${q) cc \F[A]{q)\ 



cl 
2 



(11) 
(12) 



The one-photon counting rate R v cl J corresponds to the 
classical Fraunhofer diffraction pattern. In addition, 
when q = q' , the two-photon counting rate is simply the 
square of the one-photon counting rate, 



R^(q,q)= R { *\q) 



(13) 



(6) 



On the other hand, in the case of G(x — x') = S(x — 
x'), a pair of signal and idler photons passes together 
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through the same point x at the object, and Eq. {7J) can 
be reduced as follows: 



F(q,q') = J dxA(x) 2 exp[i(q + q')x] 
1 --F[A 2 ](q + <fl, 



2tt 



(14) 



In this case, the two-photon amplitude is not separa- 
ble into any product of two independent amplitudes. In 
other words, it is a spatially entangled state. Substitut- 
ing Eq. O into Eq. 0, we get the biphoton counting 
rate: 



R^(q,q')o:\^[A 2 }(q + q')\ 



(15) 



Although the previous report on two-photon diffraction 
|2(I |2l[ discussed the case in which photon pairs are de- 
tected at separate points (q ^ q f ), we focus on the case 
in which photon pairs are detected at the same point 
(q = q') in terms of high-resolution imaging beyond the 
classical diffraction limit. Then, the two-photon counting 
rate can be rewritten as 



fl< 2 > (<?,<?)« \f[A 2 ](2q)f 



(16) 



where 2q indicates the effective transverse wave number 
of the biphoton amplitude. This result indicates that a 
photon pair behaves as a single particle with the trans- 
verse wave number 2q. On the other hand, substituting 
Eq. {E3> into Eq. (0, the one-photon counting rate be- 
comes constant: 



R {1) (q) = const. 



(17) 



This result means that the one-photon diffraction- 
interference of the biphoton will exhibit no modulation. 

From Eqs. ((12)) and (|T3|) . we understand quite natu- 
rally that these results in the case of G(x — x') = const, 
are compatible with classical optics. However, compar- 
ing Eqs. with jEJ and JE) with JEJ), we see that 
the results for the biphoton amplitude are quite differ- 
ent from those of two independent photon amplitudes in 
both one- and two-photon counting rates. Biphotons will 
exhibit half the modulation period that they would have 
in the classical case in the diffraction-interference pat- 
tern for two-photon detection. Moreover, in one-photon 
detection, no intensity modulation will be exhibited. 

To illustrate the effect of the spatial correlation on the 
one- and two-photon diffraction-interference, we calcu- 
lated and for some specific examples. In the 
following, we assume the rectangular transmission pro- 
file of the grating as shown in Figure ^ Note that in this 
case the Fourier transform of A is 



T[A\{q) 



A sm.{Ndq/2) sin(s<?/2) 



2tt sin(dq/2) q/2 



(18) 



First, we start with the case of G(x — x') = const., 
that is, an uncorrelated case. Substituting Eq. (|TH|) into 



AC*) 
A 



FIG. 1: (Color online) Assumed transmission amplitude 
through the grating, i.e., the grating function, as a function of 
the transverse position x. Ao is the amplitude transmittance 
through the grating, and d and s are the slit period and the 
slit width, respectively. The number of slits was assumed to 
be N. 
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FIG. 2: (Color online) (a) Density plot of R {2 \q,q) in the 
case of G(x — x') — const. In this plot, the two-photon count- 
ing rate R^ 2 \q^q) is expressed by a gray scale: the brighter 
color indicates larger . (b) (Lower graph) Calculated two- 
photon counting rate R^ along q = q' . (Upper graph) Cor- 
responding one-photon counting rate R^ . 



Eq. dJ, the two-photon counting rate of two indepen- 
dent photons can be obtained as follows: 



R {2) (^ 



s'm(Ndq/2) sin(sq/2) 



sin(dg/2) q/2 
sm(Ndq'/2) sin(sg'/2) " 



sin(dg'/2) q'/2 



(19) 



Figure Ufa) shows the density plot of R^ 2 \q^q') in the 
case of G(x — x') = const. The diffraction transverse 
wave numbers are normalized by q = 2ir/d, and the ratio 
d/s is assumed to be 3.2, where d and s are the slit period 
and the width of the grating, respectively. In this figure, 
the positions satisfying q = q' are indicated by a solid 
line. The lower graph in Figure |2fb) denotes the cross- 
sectional plot of R^ along q = q' '. The upper graph 
in Figure (2{b) represents the one-photon counting rate 
R {1) (q) obtained from RW(q,q') and Eq. ®. In this 
case, as expected from Eqs. (|l()j) - 0, the interference 
patterns are the same as those obtained by classical op- 
tics. 

Next we consider the case of finite spatial correlation. 
In general, it is difficult to calculate Eq. Q analytically. 



4 




FIG. 3: (Color online) (a) Density plot of R {2 \q,q) assum- 
ing a Gaussian of G(x — x). (b) (Lower graph) Calculated 
two-photon counting rate R^ along q = q' . (Upper graph) 
Corresponding one-photon counting rate 
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FIG. 4: (Color online) (a) Density plot of R {2) (q,q) in the 
case of G(x — x') = S(x — x'). (b) (Lower graph) Calculated 
two-photon counting rate along q — q' . (Upper graph) 
Corresponding one-photon counting rate 



Therefore, we may need to carry out numerical calcula- 
tions. Here we assume a Gaussian function as the trans- 
verse correlation function, 



G(x — x') 



exp 



(2V\n2)- 1 rd 



(20) 



where r is a dimensionless parameter that expresses the 
extent of the transverse correlation between the two 
photons, and corresponds to full width at half maxi- 
mum (FWHM) of the assumed Gaussian. Figure Efa) 
shows the density plot of R^ 2 \q^q') assuming r = 1. 
Comparing Figure Efa) with Figure Efa), we understand 
that the introduction of the transverse spatial correlation 
G(x — x') results in the convolution of the interference 
peaks with g(q — q') : where g(q — q') corresponds to the 
Fourier transform of G(x — x'). As a result, the peaks 
of R^ that are not present in Figure E{b) appear at 
q = q' = ±0.5, ±1.5, ... , as shown in the lower graph of 
Figure Efb). In addition, intensity modulations in R^ 
disappear as a result. Figure Ufa) presents R^ 2 \q^q') in 
the case of G(x — x') = S(x — x'), which corresponds 
to perfect spatial correlation. In this case, interference 
peaks become constant. 

Here we show the numerical calculation results for var- 




- R {2 \q,q) 1 


ur 01 




A 0.02 


i 










1.0 


. .A.I 


2.0 

. A. . 



-2-10 1 2 

q 

(b) 



FIG. 5: (Color online) Diffraction-interference patterns with 
variation of the ratio of distribution of the spatial correlation 
to the slit period (r) in (a) one-photon counting rate R^\q) 
and (b) two-photon counting rate R (y2 \q^q). 



ious values of spatial correlation width r. Figure [HI shows 
the variation of diffraction-interference patterns in both 
one- and two-photon detection. From these graphs we see 
that the diffraction-interference patterns in the region in 
which r is smaller than 0.1 is identical to that in the case 
of the perfect correlation (G(x — x') = S(x — x')). On 
the other hand, the diffraction-interference patterns in 
the region in which r is grater than 2 becomes closer to 
that in the classical case (G(x — x') = const.). Namely, 
the relationship between the extent of the transverse cor- 
relation and the object size plays a crucial role in forming 
the diffraction-interference patterns of the biphotons. 



III. EXPERIMENT 

Figure shows the schematic drawing of the exper- 
imental setup. Pairs of frequency-degenerate photons 
were generated collinearly by spontaneous parametric 
down-conversion (SPDC) in a 5-mm-long KNbC>3 crys- 
tal pumped by the second harmonic light (50 mW) of a 
single longitudinal mode Ti:sapphire laser operating at 
861.5 nm. The photon pairs were diffracted by a trans- 
mission grating (slit width: 125 /im, period: 250 /im). A 
pair of signal and idler photons has an uncertain emis- 
sion angle, which originates mainly from the divergence 
of a pump beam. Therefore, signal and idler photons may 
pass through the grating slit separately. In order to avoid 
this case, we placed the transmission grating just after 
the KNbOs crystal as shown in FigureEJb). In this geom- 
etry, each photon pair passes together through one of the 
grating slits [l^ . We separated the photon pairs from the 
pump beam by using a polarizing beam splitter (PBS) 
and a long-pass filter (LPF). By rotating a mirror (Ml), 
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FIG. 6: (Color online) (Upper figure) Schematic experimental 
setup to observe biphoton diffraction and interference. PBS: 
polarizing beam splitter, LPF: long-pass filter, Ml-2: mir- 
rors, BS: non-polarizing beam splitter, IF: interference niters, 
APD1-2: avalanche photodiodes. (Lower figure) Schematic 
drawing of the light source used for the biphoton diffraction- 
interference experiment with (b) perfect and (c) finite spa- 
tial correlation. In the case of perfect spatial correlation, we 
placed the transmission grating just after the KNbOs crystal. 
On the other hand, in the case of finite spatial correlation, we 
put the two lenses (focal length; f = 200 mm) before and after 
the KNbOs crystal in order to increase the spatial correlation 
width at the grating. 



we recorded spatial diffraction-interference patterns by 
using a two-photon detector, which consisted of a SO- 
SO non-polarizing beam splitter (BS) and two avalanche 
photodiodes (APD; EG&G SPCM-AQ161) followed by a 
multiple-stop time analyzer (EG&G 9308). To measure 
two-photon coincidence counting rates, we recorded the 
number of start-stop events within the time window (1 
ns). In addition, we simultaneously recorded the number 
of start pulses as a one-photon counting rate. In front 
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FIG. 7: (Color online) Diffraction-interference patterns of 
the parametric down- converted photons observed by (a) one- 
photon detection and (b) two-photon detection. 
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FIG. 8: (Color online) Diffraction-interference patterns of 
Tksapphire laser observed by (a) one-photon detection and 
(b) two-photon detection. 



of each APD we placed an interference filter (IF; cen- 
ter wavelength A c = 860 nm, bandwidth AA = 10 nm). 
To compare the results with those of classical light, we 
also observed, using the same apparatus, the diffraction- 
interference patterns of the Ti:sapphire laser and that of 
thermal light generated from a tungsten-halogen lamp. 

Finally, we observed the diffraction-interference pat- 
terns in the case of finite spatial correlation. To achieve 
this, we observed diffraction-interference by inserting two 
lenses (focal length; / = 200 mm) before and after the 
KNb03 crystal and placing the grating at a distance of 
40 cm from the crystal as shown in Figure EJc). The 
increased divergence of the pump beam enlarges the un- 
certainty of the transverse wave number of the SPDC 
photons. This results in the increase of the spatial corre- 
lation width between the two photons at the grating . 
Therefore, by using the lenses we can control the spatial 
correlation between the signal and idler photons. 
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FIG. 9: (Color online) Diffraction-interference patterns of the 
halogen lamp observed by (a) one-photon detection and (b) 
two-photon detection. 
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FIG. 10: (Color online) Diffraction-interference patterns of 
the parametric down- converted photons observed by (a) one- 
photon detection and (b) two-photon detection when two 
lenses are inserted before and after the KNbOs crystal. 



IV. RESULTS AND DISCUSSION 

Figures and |HI show the measured diffraction- 
interference patterns of parametric down- converted pho- 
tons and those of the Tiisapphire laser, respectively. In 
both figures, the upper graphs represent the interference 
patterns observed by one-photon detection and the lower 
graphs are those observed by two-photon detection. The 
open circles in each figure represent the measured data 
points. The transverse wave numbers are normalized by 
qo = 27r/cZ, where d is the slit period of the grating. These 
experimental data were fitted with the theoretically ex- 
pected functions (fT2 |) -(|13 |) . (|lfi )) -(|17 )) . as indicated by the 
solid curves. One can see that the experimental data are 
in close agreement with the theoretical prediction |22| • In 
particular, the two-photon interference of the SPDC seen 
in FigureEIb) exhibits half the modulation period of that 
of the Ti:sapphire laser in Figure |BI (see Eq. (|16[)). This 
result indicates that the biphoton generated by SPDC be- 
haves as an effective single particle that is associated with 
half the wavelength of the constituent photons. In other 
words, the diffraction of a biphoton can be explained by 
the concept of the photonic de Broglie wave. 

It is also noteworthy that the one-photon interference 
of the SPDC exhibits no modulation, whereas that of the 
Ti: sapphire laser exhibits normal modulation that can be 
understood by means of classical optics. In classical op- 
tics, the disappearance of an interference fringe might be 
caused by the shortage of coherence length owing to the 
wide spectral distribution of the parametric emission. To 
be sure that the observed phenomenon originates from a 
quantum mechanical effect but not from classical effects 
caused by the spectral width of the light, we also mea- 
sured the diffraction-interference pattern of thermal light 
from a halogen lamp, as shown in Figure 03 Here the 
thermal light was put through a pinhole with a diameter 
of 50 fim to improve spatial coherence. We note that 
the spectrum of the SPDC emission at the detectors is 
almost the same as that of the halogen lamp, because we 
detected both emissions through the same interference 



filters. However, the measured pattern of the halogen 
lamp is quite different from that of the SPDC, and is 
quite similar to that of the Ti:sapphire laser. Thus the 
disappearance of the one-photon diffraction-interference 
pattern in the SPDC emission is not due to its spectral 
width. The disappearance of the one-photon interfer- 
ence pattern is explained by the lack of spatial coher- 
ence, which inevitably originates from the spatial corre- 
lation of the photons. To discuss the spatial coherence of 
the constituent photon in a spatially correlated state, we 
consider the first-order coherence function defined by 



g {1) (x,x') = {^\d\x)d(x')\^) . 



(21) 



We assume that the input state is the two-photon 
state given in Eq fQ, and the two-photon probability 
amplitude F(x,x f ) is given by the following perfectly cor- 
related state: 



F(x,x')=A(x)A(x')6(x-x'). 



(22) 



Substituting Eqs. JIJ and into Eq. l|2T|l. we get the 

first-order coherence function: 

g (1) (x,x') = 2 (\A(x)\ 2 + \A(x')\ 2 ^j A*(x)A(x')S(x - x 1 ). 

(23) 

This result indicates the constituent photons in the spa- 
tially correlated state are completely incoherent. As a re- 
sult, the one-photon interference pattern is washed-out. 
The wide distribution of each transverse wave number 
q (or q') as shown in Figure Q] is another aspect of the 
spatial incoherence. Thus the disappearance of the one- 
photon interference pattern is a natural consequence of 
the spatial correlation. 

Next we show the observed diffraction-interference pat- 
terns in the case in which the two lenses are placed be- 
fore and after the KNbOs crystal, as shown in Figure ITU1 
Comparing these results with the former ones shown in 
Figures. [7| and [HI we see that the diffraction-interference 
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pattern in the two-photon detection method (Fig. lTDT b)) 
exhibits the intermediate pattern between the two ex- 
treme cases. Furthermore, the biphotons begin to show 
partial modulation in the one-photon detection method 
(Fig. irUT a)). again exhibiting the intermediate pattern 
between Figures [3 and [HI Considering a finite correla- 
tion width, we can also reproduce these patterns from 
Eqs. Q)-©- Here we assume the Gaussian correlation 
function: 



G(x — x') — exp 



x — x' V 
0.56d ) 



(24) 



which corresponds to r = 0.78. The calculated pat- 
terns are in close agreement with the measured patterns 
as shown by the solid curves in Figure EH Although 
some previous experimental reports have focused on spa- 
tial correlation in for ming the diffraction patterns of a 
Young's double slit p3, 24], we have described the present 
experiments in terms of high-resolution imaging using 
spatially correlated photons and also shown the validity 
of the Fourier-optical analysis taking account of the spa- 
tial correlation between constituent photons. In terms 
of such high-resolution imaging as quantum lithography, 
these results indicate that the distribution of a transverse 
correlation has to be much smaller than an object so as 
to achieve a spatial resolution high enough to surpass 
the classical diffraction limit. Otherwise, even entangled 
photons behave as a classical light source according to the 
relationship between the extent of the transverse correla- 
tion and an object. However, there is little advantage in 
using SPDC as an entangled photon source for quantum 
lithography because the spatial resolution of a pair of 
parametric down- converted photons is the same as that 
of the pump beam for SPDC. One of the methods to over- 
come this problem is to use Hyper Parametric Scattering 
(HPS) as discussed by Strekalov |25| . The spatial reso- 
lution of entangled biphotons generated by HPS should 
be greater than that of the pump beam by a factor of 
2, whereas the wavelength of the generated photons is 
equal to that of the pump beam. By utilizing HPS, in 
the future it will be possible to obtain a spatial resolu- 
tion of 100 nm or below by using ultraviolet laser light. 
Recently we have successfully generated ultraviolet en- 
tangled photons, whose wavelength is approximately 390 
nm, via HPS in a semiconductor crystal [26|. This may 
open the way to quantum imaging with truly high reso- 
lution beyond the classical limit. 

Our Fourier-optical analysis is applicable to various 
types of two-photon diffraction experiment, such as 
"ghost interference" |5|. The conditions of the ghost in- 
terference experiment differ from those of our experiment 
in regard to two points. First, the constituent photons 
of a photon pair are divided into two arms: the photon 
amplitude at q (or q') consists of only one photon, whose 
position is denoted by x (#'), scattered by the object A 
(B). Second, the position of one detector placed after the 
object is fixed, while the other is scanned. Under these 
conditions, the term of A(x')B{x) in Eq. (0) contributes 



zero, and thus the two-photon amplitude at the object 
becomes 



F(x,x') = -A(x)B(x , )G(x ■ 



(25) 



In addition, the creation and annihilation operators at 
the object are commutable because the two photons do 
not exist on the same object. As a consequence, the two- 
photon amplitude at the far field is given by 



Hq>J) = ^ J dx J dx'A(x)B(x')G(, 



x — x') 



(26) 



x exp [i(qx + qx)] . 



Since the object is placed only in one arm, while there is 
no object in the other arm, 



1 

47T 



dx / dx' A(x) 



x G(x - x) ex.p[i(qx + qx)}. (27) 



The diffraction pattern is given by R i2) (0,q f ) in Eq. ®, 
which corresponds to the cross-sectional plots along q = 
0. In the classical case of G(x — x') = const., Eq. (|2?|) 
becomes 

F(q,q') = ^- J dx J dx'A(x) ex.p[i(qx + q'x')) 

= - J dxA(x)exp(iqx)5(q'). (28) 



Therefore, 



i? (2 W) oc \T[A\{q)5{ q ') 



oc 



R$(q) l%')l 2 . 



(29) 



Note again that J 7 [A] represents the Fourier transform of 
A. Thus there is no diffraction pattern as a function of 
q' (Fig. HIT a)). On the other hand, in the quantum case 
of G(x — x') = S(x — x'), Eq. (|27jl reduces to 



i 

47T 

oc T[A](q + q'). 



Thus, in accordance with Eqs. (g 
photon counting rate is given by 



dxA(x) exp[z(<2 + q')x\ 

(30) 

| and O, the two- 



fl< 2 W) oc \T[A]{ q + q')\ 2 



(31) 



This result indicates that the diffraction pattern pre- 
sented by spatial entangled photon pairs appears as a 
function of q' . As explained in the preceding section, we 
can also understand this image formation process as the 
diffraction pattern spread in the direction perpendicular 
to the line satisfying q = q' as a result of the spatial 
entanglement (Fig. HTT b)). 
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FIG. 11: (Color online) Density plot of R^ 2 \q,q) in the case 
of the ghost interference experiment, assuming a grating as 
the object, (a) Classical case, that is, G(x — x') — const, 
(b) Quantum case, that is, G(x — x') = S(x — x). (c) A 
density plot when the plot (a) oscillates in a direction parallel 
to the line satisfying q = q' . (d) A comparison of the ghost 
diffraction patterns. The upper graph is the cross-sectional 
plot along q = seen in Figure (c) and the lower one is that 
in Figure (b). 



Recently, several groups have reported coincidence 
imaging experiments using classical light sources, which 
stimulated discussion regarding the difference between 
quantum and classical coincidence imaging j^J [53, H3, 
Efj l3(]| . Two types of experiment use classical light 
sources. One is the experiment using a classically 
momentum-correlated bea m l27l l28l. The other is that 
using thermal light sources Here we dis- 

cuss only the experiments using a classically momentum- 
correlated beam in terms of our Fourier-optical analy- 
sis because a discussion taking into account a statisti- 
cal average is required for incoherent illumination, such 
as in the experiment using thermal light. The most 
significant feature of the experiment using momentum- 
correlated beams is that the diffraction pattern of clas- 
sically momentum-correlated beams, produced by rotat- 
ing a mirror before the beam splitter as shown in Fig- 
ure [EH a )? appears as a function of q' whereas there is no 
spatial entanglement on the object, that is, G(x — x') = 
const. We can understand this image formation process 
as follows. As shown in Figure ITTT a), the intensity mod- 
ulation of classical beams appears as a function of q but 
not q' . However, by rotating the mirror, the diffraction 
pattern shown in Figure ITTT a) oscillates in a direction 
parallel to the line satisfying q = q' with the scanning 
frequency of the mirror. Therefore, taking an average of 



grating 




FIG. 12: (Color online) (a) Schematic view of the ghost in- 
terference experiment performed in Ref. [2^. Here we put a 
grating as the object instead of a double slit, (b) An assumed 
classical version of the experiment described in this paper. 
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(a) (b) 

FIG. 13: (Color online) (a) A density plot when the R {2) (q, q) 
in Figure Ufa) oscillates in a direction parallel to the line sat- 
isfying q — q[. (b) (Upper) A cross-sectional plot along q — q[ 
and (Lower) that along q — —q' in the plot (a). 

the product of the two signals over a sufficiently long pe- 
riod compared to the scanning period of the mirror, we 
can obtain the two-photon counting rate as follows: 

# (2 W) = [ R ( k\q-qo)R ( i ) (q' -qo)dq , (32) 

where qo represents the displacement of the pattern by 
the rotation of the mirror, and and R^ are the one- 
photon counting rate from each object having the form of 
A(x) and B(x), respectively. Since R^\q f — qo) is given 
by the delta function 5(q' — qo) as a result of B(x') = 
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const, in the ghost diffraction experiment, Eq. can 
be reduced as follows: 

R (2) {q,q') = J R ( l\q ~ ?o)%' - <?o)%>, 

= R^(q-q'). (33) 

The density plot of Eq. is shown in Figure ITTT c). 
This is a symmetric plot of that in Figure ITTT b) with 
respect to q = 0. Thus, the cross-sectional plots along 
q = seen in Figures ITTT b) and (c) give the same pat- 
tern, as shown in Figure lTTT d). Therefore, to identify the 
difference between the pattern produced by quantum cor- 
related photons and that by classically correlated beams 
in the case of ghost interference experiments, more de- 
tailed discussion is required as reported in Ref. HH • 
On the other hand, the difference clearly appears in the 
case of the experiment in which two photons pass through 
the same object, that is, R^\q) = R^\q) in Eq. (|32|) . 
An assumed schematic setup of the classical version of 
this kind of experiment is shown in Figure H^T b). when 
the form of the transmission amplitude A{x) is the same 
as that in Figures EMI Figure ITST a) shows the expected 
density plot, which corresponds to the pattern shown in 
Figure EJa), a pattern that oscillates in a direction par- 
allel to the line satisfying q = q' '. In addition, the upper 
and lower graphs in Figure ITST b) represent their cross- 
sectional plots along q = q' and q = — respectively. 
Interference peaks appear with half a modulation period 
of the normal diffraction-interference pattern shown in 
Figure El in the case of q = q', while it shows no modula- 
tions in the case of q = q' . To clarify the image formation 
process, we rewrite Eq. as follows: 

R {2) (A) = J rV> (A + Q)R<U (Q)dQ, (34) 

where 

Q = q'-q , A = q-q'. (35) 

If we detected at the same point, i.e., q = q', the two- 
photon counting rate R( 2 \A) becomes constant as a re- 
sult of A = 0. On the other hand, in the case of detection 
at the symmetrical position q = — q f , Eq. (|34|) becomes 

R?\2q) = J R { l\Q + 2q)R { l\Q)dQ. (36) 

This equation means that the interference of the two 
photons has half a modulation period of the one-photon 
counting rate R^\q). However, the two-photon count- 
ing rate R^ is nothing more than the autocorrelation 



function of the one-photon counting rate R A . Thus, the 
diffraction width of this pattern is no narrower than that 
of the normal patterns as shown in Figure EJb). This 
result is consistent with the consideration that spatial 
resolution beyond the diffraction limit can be achieved 
only by using spatially entangled photons but cannot be 
achieved using a classical light source, as discussed in 
Ref. |U. 

V. CONCLUSIONS 

The experiment described in this paper focuses on the 
use of entangled photons for novel quantum imaging tech- 
nologies, which has been attracting much attention in 
connection with the concept of the photonic de Broglie 
wave. We have investigated the role of spatial corre- 
lation in diffraction-interference phenomena by measur- 
ing the spatial diffraction-interference patterns of spon- 
taneous parametric down-converted biphotons through 
a transmission grating. We have shown that the spa- 
tially correlated biphoton exhibits half the modulation 
period of that of classical light. Also, we have confirmed 
that the one-photon interference of the biphoton exhibits 
no modulation. Furthermore, by controlling the spatial 
correlation between the two photons, we have success- 
fully demonstrated that the diffraction-interference pat- 
tern changes from the perfectly correlated biphoton case 
to the classical case. In order to understand these ex- 
perimental results in a quantitative and straightforward 
manner, we developed the Fourier-optical analysis of a 
two-photon state taking into account the spatial correla- 
tion between the two photons. 

In our experiment, we have investigated the diffraction- 
interference patterns formed by a one-dimensional reg- 
ular grating, and demonstrated the validity of Fourier- 
optical analysis. Nevertheless, we would like to empha- 
size that this analysis is extendable to arbitrary one- and 
two-dimensional objects, and can thus be is very useful 
for future imaging technologies utilizing entangled pho- 
tons. 
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